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>^ ■ In this paper we investigate regular patterns of matrix elements of the nuclear shell model Hamil- 

tonian H, by sorting the diagonal matrix elements from the smaller to larger values. By using simple 
plots of non-zero matrix elements and lowest eigenvalues of artificially constructed "sub-matrices" 
h of H , we propose a new and simple formula which predicts the lowest eigenvalue with remarkable 
precisions. 
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The diagonalization of matrices is a fundamental practice in nuclear structure physics as well as many other fields. 
However, diagonalization becomes difficult if the dimension of the matrix is very large. Statistical approaches are 
very suggestive and have been developed, e.g., in Refs. [ll-Q: where the lowest eigenvalue is presented in terms of the 
energy centroid and spectral moments. 

Recently, we showed in Ref. [8| that sorting diagonal matrix elements of a given nuclear shell model Hamiltonian 
from the smaller to the larger values provides us with a new approach to evaluate the eigenvalues. By sorting the 
diagonal matrix elements, we are able to evaluate all eigenvalues based on a very strong linear correlation between the 
diagonal matrix elements and exact eigenvalues. This method was found to work very well for medium eigenvalues 
but deviates for the lowest ones. However, in nuclear structure physics as well as many other fields of sciences, we 
are interested in the low-lying states. It is therefore very desirable to refine the approach towards more and more 
accurate evaluation of the low-lying eigenvalues, by sorting the diagonal matrix elements. 

In this paper, we propose a new approach to predict the lowest eigenvalue of the nuclear shell model Hamiltonian. 
To exemplify our method, we shall use a few realistic examples of nuclear shell model calculations. All results in this 
article are based on the shell model code by the Kyushu group @t11| ■ The shell model basis states of this code are 
constructed by using the coefficients of fractional parentage discussed in Ref. [ll[. In this paper we take the USD 
interactions of Ref. [12] ■ Other interactions such as the Yukawa-type interactions of Refs. [ifljIlSl S^'^^ similar results. 

Let us denote the matrix of spin / states of the nuclear shell model Hamiltonian H by H^^\ and the matrix 
elements of H'^^'^ by Hi- , where i and j represent indices of basis states. In Fig. [TJ we present two typical examples 

of distributions of the magnitude of H\- , based on the J'^ = 0+ and J'^ — 2+ states of the ^^Mg nucleus. The color 
from blue to red corresponds to values from zero to large magnitudes. From panels (a-b) of Fig. [U one sees that the 
values of Hi- ' (panels (a-b) of Fig. [1]) look "random" . However, if one sorts the diagonal matrix elements from the 

smaller to larger values, as in Refs. Q, the values of Hi- ' decrease rapidly and become zero if they are "far" enough 
from the diagonal line, as shown in panels (a'-b') of Fig. [T] 

Let us investigate this behavior in another form. We study the probability for H^, to be non-zero (after sorting 
the diagonal matrix elements of H^- ), while moving away from the diagonal line, versus d, denoted by p{d) = 

j)_^''+'' , d = 1, 2, 3, • • • ,D. Here d is the "distance" of H-,- ' from the diagonal line, and D is the dimension of 

matrix H^, ' for spin / states. As shown in Fig. ^a.) and (b) for the I'" = 0+ and 2+ states of the ^^Mg nucleus, p{d) 
becomes zero at a critical value d = do! the value of In dp equals 6.75 and 8.22, respectively. 

An argument for the regular patterns described in Fig. 1 and Fig. 2 (a) and (b) is as follows. With the diagonal 
matrix elements sorted from the smaller to larger values, one classifies configurations from the lowest to the largest in 
energy, roughly by particle-hole excitations. Configurations that come first are the lowest, and the states that come 
last are n-particle-n-hole excitations out of those low configurations. Because the shell-model Hamiltonian consists of 
one-body and two-body operators, one can not connect those configurations that are "distant" , e.g., the configurations 
with the lowest energy and n-particle-n-hole configurations with n > 2. This explains the reason why all values of 
i?j- become zero for d > do. Soon we shall find that the value of do is very important in predicting the lowest 

eigenvalue of the matrix if '^^^. 

In Refs. [3,0,01 the lowest eigenvalue is presented in terms of InD, where D is the dimension of the matrix ij'^^ 
Although the formulas of the lowest eigenvalues presented in Refs. 0, |a, 0] are applicable to the random ensemble 
average (not to individual sets of interactions parameters) , one naturally asks whether or not certain plots of the lowest 
eigenvalue versus the dimension could be useful in evaluating the lowest eigenvalue of realistic systems studied in this 
paper. Let us sort the diagonal matrix elements from the smaller to the larger, as in Refs. [8|. Then we truncate 
artificially the matrix i?'^^ and obtain a "sub-matrix" h with dimension d {d < D), and hij = ij|- ' (i, j = 1, 2, • • • , d). 
We diagonalize h and obtain the lowest eigenvalue e^ of the matrix h, and plot e^ versus Ind. In Fig. 2(a'-b') we 
present the e^-lnd plots for the /'^ = 0+ and 2+ states of the ^"'Mg nucleus. One sees that e^ decreases linearly with 
Ind when d is smaller than a critical dimension d = Dq, and decrease again linearly with Ind but with a smaller 
slope. Apparently, the value of Dq and the slopes for both d < Dq and d > Dq suffice for the evaluation of the lowest 
eigenvalue of Hl^ . 

In Fig. 2 one sees that the values of do where jo(do) = in Fig. 2(a) and (b) coincide with Dq in Fig. 2(a') and 
(b'), respectively. Panels (a,b) are based on the same matrices as (a',b'), respectively. For convenience, we use the 
same scale in panels (a, a') and (b,b'), and plot two dotted lines to guide the eyes in order to see such coincidence. 
From Fig. 2 one also sees that the slope for d > do (denoted by k') is smaller than that for d < do (denoted by k). 

An intuitive understanding of the facts that do — Dq and fc' < fc is given as follows. Because H^- ' are zero when 
d > do (i.e., p{do) = for i,j > do), there is no contribution to the lowest eigenvalue from these matrix elements. On 
the other hand, some of matrix elements i?! and H^^ , with < i < D — do and j > D — do, are non-zero (see Fig. 
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FIG. 1: (Color online) Magnitude of the matrix elements without (panels (a) and (b)) and with (panels (a') and (b')) sorting 
the diagonal matrix elements from the smaller to the larger. The color from blue to red corresponds to values from zero to large 
magnitudes. The results are based on I^ — 0"*" (panels a and a') and I^ = 2"*" (panels b and b') states of the Mg nucleus, 
obtained by using the USD interactions. The magnitude of Hi- without sorting the diagonal matrix elements (left hand side) 
are close to "random" , and those with sorting the diagonal matrix elements (right hand side) decrease rapidly as going farther 
from the diagonal line. 

l(a') and (b')), lowering down the smallest eigenvalue of the h matrix. The slope of the erf-lnd plot therefore changes 
at d = do and becomes smaller for d > do than that of h ior d < dg. 

In Fig. 3 we present our results of \n{Do) /\n{do) and fc'/fc, based on 200 examples of ij'^^ for a number of sd 
shell nuclei by using the USD interactions. One sees that ln(Do)/ln(do) are very close to 1.0 and that k' ~ |A: with 
fluctuations. For simplicity we assume that ln(Do)/ln((io)=l and k' = -rk for all cases throughout this paper. 

Making use of these regularities, we obtain a new and simple formula to evaluate the lowest eigenvalue of the matrix 
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where k and b are the slope and intercept of the e^-hid plot for d < do; do is determined by p{do) = 0, and D is the 
number of spin / states. Because the e^-lnd plot shows a nice linearity (see Fig. 2(a',b')), we extract the values of k 
and b based on sub- matrices h of H'^^^ with d < D/IO for all cases. 

In Fig. 4 we present a comparison of the lowest states of spin I predicted by the above Eq. ([1]), and those obtained 
by the linear correlation (i.e., Eq. (3) of Ref. Q), with those calculated by diagonalizing ij'^^ (/ = 0,2,4,6,8) for 
two nuclei, ^''Mg and ^®Si. One sees the remarkable agreement between the exact eigenvalues (the column "exact") 
and our predicted ones ("predl"), and substantial improvements achieved by Eq. ([1]) in comparison with Eq. (3) of 
Ref. Q ("pred2" in Fig. 4). Without going into details, we mention that the overall root- mean-squared deviation £ 
for the two-hundred cases we checked in Fig. 3 (defined by £" = ^^ ^^^cxact _ ^prcd^2^^^ where N is the number 
of examples that we checked and here N = 200) is 0.38 MeV, assuming that do = Do, k' = |A: for all examples. 

Here we give a brief discussion of formulas in Refs. 0, |j, M^ and the formula proposed in this paper. Ref. [J] 
reported a correlation between the lowest eigenvalue and the spectral width a of the spin / states. Ref. [3] suggested 
a simple formula E^(^ = Ej — ^ln£)/ln2 a, where Ej is the average energy of spin / states. Ref. [y] suggested a 
formula £'^^„ = Ej - y/alnD + b a (similar to Ref. Q), with a ~ 1.00 and b ~ 0.40. Ref. H] refined the results 
of Ref. [y] by including the third moment analytically, with an additional factor (1 — -^-^ (^) 1 multiplied in the 




FIG. 2: (Color online) Distribution for Hl^ to be non-zero, denoted by p{d), and the lowest eigenvalues of h, denoted by 
td{h). The results of (a,a') and (b,b') are obtained for the J^ = 0^ and J"^ = 2"*" states of ^^Mg, respectively. One sees that 
the values of d where p(do) = and where the td{h)-\nd plots change the slope (-Do) coincide, i.e., do = Dq. We use the same 
scale in (a,a') or (b,b'). The dotted lines are used to guide the eyes. 



second term of the formula in Ref. (6[. Here (T3 is the third moment of the eigenvalues. These formulas arc applicable 
to the random Hamiltonians statistically (i.e., the ensemble average), not to the individual Hamiltonians such as 
for realistic systems. In Refs. [8|, the formula by using the linear correlation between the eigenvalues and diagonal 
matrix elements is applicable to individual sets of parameters and works well for medium eigenvalues, but it docs not 
work very well for the lowest (or the largest) eigenvalues. The formula proposed in this paper is found to predict 
remarkably well the lowest eigenvalues of the nuclear shell model Hamiltonian, as shown in Fig. 4. We should also 
add that there are many other efforts towards overcoming the limitation of dimension in diagonalizing large matrices, 
see refs. [li-[l3- 

To summarize, in this paper we first investigated the distribution of non-zero off-diagonal matrix elements of the 
nuclear shell model Hamiltonian. We demonstrated that the non-zero off-diagonal matrix elements exhibit regular 
patterns, if one sorts the diagonal matrix elements from the smaller to larger values; without sorting the diagonal 
matrix elements, the off-diagonal matrix elements look random. Almost all matrix elements becomes zero, if the 
matrix elements are "distant" enough from the diagonal line, after sorting the diagonal matrix elements. 

A very simple formula of the lowest eigenvalue for the shell model Hamiltonian matrix H^-^' is proposed, based on 
the regular patterns of the p{d) and ed-lnd plots for sub-matrices h. There exists a "critical" dimension, Dq, at which 
the slope of the Cd-lnd plot changes. The slope for d> Dq is empirically found to be equal to 2/5 of that for d < Dq 
with fluctuations. The value of Dq is found to be equal to do which can be obtained easily by using p{do) ~ 0. Here 
p{d) represents the probability for Hij to be non-zero while moving away from the diagonal line. 

The overall root-mean-squared deviation for two-hundred shell model Hamiltonians of nuclei in the sd shell is 0.38 
MeV (the relative deviation is about 0.003), assuming that do = Dq and k' = ^k, with k and b obtained from the 
sub-matrices h of H^^K The dimension of h is much smaller than D, and here we take d < D/10. This demonstrates 
that our predicted results of the lowest eigenvalue based on our new formula are in very good agreement with the 
exact values, even if one treats much smaller "sub-matrices" of H^^\ We therefore expect that our new formula has 
significance for future theoretical studies of nuclear structure. It will be also interesting to investigate whether or not 
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FIG. 3: (Color online) Distribution for ln(_Do)/ln(do) and k'/k for two hundred J/'^^ matrices for a number of nuclei in the 
sd shell, the dimension D of which goes from 500 to 10000, by taking the USD interactions. One sees that ln(Z)o)/ln(do) — 1, 
and that k' ~ |fc with considerable fluctuations. In this paper we assume that ln(Do)/ln(do) = l and k' — |fc when we predict 
the lowest eigenvalue of H^^' , for simplicity. 

other low-lying states have similar features. 

What has not been yet understood at a microscopic level is why the Cd-ln d plot exhibits a remarkable linearity. 
Further consideration of these issues is warranted in future studies. 
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FIG. 4: (Color online) Comparison of low- lying levels obtained by exact diagonalization (denoted by "exact") and those 
predicted by using Eq. ([T| (denoted by "predl") and the linear correlation formula suggested in Ref. [31 (denoted by "pred2") 
for the ^''Mg nucleus in panel (a), and the '^^Si nucleus in panel (b). The results are obtained by using the Yukawa-type 
interaction of Refs. [l0,[l3|. One sees substantial improvement of predictions for the lowest eigenvalues for each spin / by using 
Eq. (1) of this paper. 
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